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BACKGROUND 


Modem  explosively  formed  projectile  warheads,  such  as  SADARM,  produce  a  single  robust  penetrating 
fragment.  However,  greater  lethality  could  often  he  achieved  by  either  a  number  of  smaller  fragments  or  a 
deep  penetrating  jet  depending  on  the  target  properties.  Therefore,  multiple  mode  warheads  provide  increased 
utility  and  overall  lethality  due  to  their  applicability  over  a  larger  range  of  targets.  The  Advanced  Warhead 
Concepts  ILIR  project  is  investigating  advanced  multiple  mode  warhead  concepts  based  on  charge  designs 
employing  explosives  with  different  properties  and  initiation  logic.  These  advanced  warhead  concepts  are 
being  investigated  using  dynamic  finite  clement  analysis  and  experimental  verification.  Current  dynamic  finite 
element  programs,  such  as  DYNA3D,  do  not  adequately  treat  detonation  interaction  phenomena  produced  by 
explosives  with  different  properties  and  multiple  initiation  points.  Therefore,  the  FY90  effort  was  to  develop 
and  implement  a  significant  enhancement  for  the  treatment  of  detonation  interactions.  This  has  been  done  by 
producing  an  advanced  detonation  products  equation  of  state,  appropriate  for  overdriven  detonations  produced 
by  detonation  interactions.  This  new  equation  of  state  provides  an  improved  high  pressure  description,  and 
retains  the  low  pressure  expansion  behavior  required  for  standard  material  acceleration  modeling. 

APPROACH 

Thcrmochcmical  calculations  have  proven  to  be  very  useful  for  the  prediction  of  explosive  products 
properties,  particularly  near  and  above  the  Chapman- Jouguct  state.  Unfortunately,  they  do  not  reproduce  the 
products  expansion  behavior  accurately  enough  for  typical  warheads  design.  Currently,  thermodynamic 
equations  of  state  (JWL)  used  for  warheads  design  are  normally  calibrated  to  give  agreement  with  copper 
cylinder  explosive  expansion  experiments.  'Ihesc  equations  of  slate  have  not  been  calibrated  for  high  pressures 
above  the  Chapman-Jouguet  state.  Experimentation  and  comparison  with  thermochemical  calculations  (figures 
1,2,  and  3)  have  demonstrated  that  a  poor  description  of  the  high  pressure  region  exists.  In  order  to  achieve  a 
suitable  equation  of  state,  an  appropriate  equation  of  state  form  (JWLB)  was  derived.  A  standard  explosive 
(octol  75/25)  was  used  for  copper  cylinder  expansion  tests  and  explosive  dent  plate  tests.  The  new  equation  of 
state  was  parameterized  for  octol  75/25,  using  both  thermochcmical  calculations  and  cylinder  test  results  in 
order  to  adequately  describe  both  the  high  pressure  region  and  lower  pressure  expansion  behavior.  After 
implementing  the  equation  of  stale  into  a  finite  element  program,  dynamic  finite  element  analysis  of  the 
experimentation  using  JWLli  and  a  standard  equation  of  state  (JWL)  with  standard  octol  75/25  parameters  was 
completed.  The  experimental  and  computational  results  were  reduced  and  compared. 

EQUATION  OP  STATU  IORMUI  A'llON 

’Ihc  equation  of  state  form  was  chosen  so  as  to  adequately  describe  the  high  pressure  regime  produced 
by  overdriven  detonation,  and  yet  retain  the  low  pressure  expansion  behavior  required  for  standard  material 
acceleration  modeling.  To  this  end,  the  derived  form  is  based  on  the  Joncs-Wilkens-Lee  (JWL)  equation  of 
stale  due  to  its  asymptotic  approach  to  an  ideal  gas  at  high  expansions.  Additional  exponential  terms  and  a 
variable  Clruneiscn  parameter  have  been  added  to  adequately  describe  the  high  pressure  region  above  the 
Chapman-Jouguet  shite.  The  resulting  equation  of  suite  form,  named  Jones-Wilkens-Lee-Bakcr  (JWLB),  is  as 
follows: 


The  JWLB  equation  of  state  form  is  based  on  a  first  order  expansion  around  the  principle  isentropc: 

Ps  -  jA.e*RiV  +  CV^1). 

Using  the  Gruneisen  Parameter, 


k  =  V— 
dE 


V 


the  isentropic  identity, 


p  — .  HE. 

dV 

and  the  Chapman-Jouguet  condition, 


Ecj  =  Eo  +  j(Pcj  +  Po)(Vo  -  Vcj) 


the  final  form  may  be  derived 


P  =  — (E  -  Es)  +  Ps  *  Final  Form. 

Some  important  characteristics  of  the  equation  of  state  are  that  the  Gruneisen  Parameter  A,  is  represented  as  an 

V  <?P 

analytic  function  of  specific  volume,  V.  A  +  1  approaches  a  constant  adiabatic  gamma, - -  to  +  1 

PiVs 

for  large  V,  so  that  ideal  gas  behavior  is  asymptotically  approached. 

EXPERIMENTATION 


Copper  cylinder  expansion  tests  and  steel  dent  plate  tests  were  done.  The  copper  cylinder  expansion 
tests  consisted  of  oxygen  free  copper  tube  (10"H  X  POD  X  3/4TD)  filled  with  octol  75/25.  Each  charge  was 
detonated  from  one  end,  and  the  cylinder  wall  movement  was  measured  using  an  ultrahigh  speed  smear 
camera.  The  dent  plate  tests  consisted  of  a  3"  diameter  charge  placed  on  top  of  three  3”  thick  RHA  steel 
plates.  The  charge  is  detonated  and  the  resulting  dent  profile  is  measured  using  a  sliding  micrometer.  The 
charges  consisted  of  1.5"  and  2.5"  height  octol  75/25. 

PARAMETERIZATION 


A  generalized  method  of  parameterization  has  been  developed  using  a  nonlinear  optimization  program 
(NLQPEB),  thermochemical  calculations  (TIGER),  and  a  mathematical  cylinder  test  model  (CYLTEST).  The 
NLQPEB  program  was  developed  by  ARDEC  for  the  solution  of  generalized  nonlinear  optimization  problems 
with  equality  and  inequality  constraints.  The  program  utilizes  sequential  quadratic  programming,  with  a 
Broyden-Fletcher-Goldfarb-Shanno  (BFGS)  variable  metric  update  and  a  modified  merit  function  based  on  the 
Kuhn-Tucker  first  order  necessary  conditions.  NLQPEB  is  used  to  solve  for  a  set  of  JWLB  parameters  that 
reproduce  the  predicted  pressure,  volume,  Gruneisen  parameter  principle  isentrope  behavior  in  the  least  squares 


sense.  Additional  equality  constraints  are  imposed  so  that  the  predicted  Chapman-Jouguet  state  be  reproduced, 
and  that  CYLTEST  reproduces  the  experimental  cylinder  expansion  test  results  at  seven  volume  expansions. 
The  pressure,  volume,  Gruneisen  parameter  principle  isentrope  behavior  and  the  Chapman-Jouguet  state  are 
predicted  using  the  computer  program  TIGER.  Either  the  Becker-Kistiakowsky-Wilson  (BKW)  or  the 
Jacobs-Cowperthwa  ite-Zwisler  (JCZ3)  thermochemical  equations  of  state  can  be  used  for  the  calculations.  The 
mathematical  cylinder  expansion  test  model,  implemented  into  the  program  CYLTEST,  is  a  modification  of  the 
model  suggested  by  G.I.  Taylor  using  Reynolds  hydraulic  treatment.  The  modification  includes  radial  flow  and 
cylinder  wall  thinning  considerations.  The  parameterization  method  was  used  with  BKW  and  the  LLNL 
revised  BKW  parameters  (BKWR)  to  produce  JWLB  parameters  for  octol  75/25.  The  octol  75/25  JWLB 
parameters  are: 
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Pcj.Mbar 

1.821 

.096034 

.85381 

.33755 

Figures  4,  5,  and  6  compare  BKWR  and 

JWLB  calculations. 

DYNAMIC  FINITE  ELEMENT  COMPUTATIONS 

The  equation  of  state  has  been  implemented  into  the  dynamic  finite  element  program  DYNA2D. 
The  cylinder  expansion  and  plate  dent  tests  have  been  modeled  using  JWLB,  as  well  as  JWL  with  standard 
parameters  for  comparison.  Figure  7  shows  the  cylinder  lest  modeling  at  6  microsecond  intervals.  Figure  8 
shows  the  plate  dent  test  initial  and  final  configurations.  The  material  boundaries  and  computational  mesh  of 
both  computations  can  be  clearly  seen  in  the  figures.  The  Johnson-Cook  material  model  and  the 
Mic-Gruneisen  equation  of  state  were  used  for  both  the  copper  and  steel  material  descriptions.  Extensive 
rezoning  was  required  for  the  plate  dent  calculation. 

RESULTS  COMPARISON 

Figure  9  presents  a  comparison  of  the  finite  clement  and  experimental  cylinder  expansion  test 
results.  The  cylinder  test  dynamic  finite  element  results  agree  very  closely  with  the  experimental  results.  The 
JWLB  equation  of  state  results  agree  slightly  better  than  the  JWL  results.  The  dent  plate  computations  also 
agree  closely  with  the  experimental  results.  It  was  found  that  the  closest  agreement  was  achieved  using  an 
Ignition  and  Growth  reactive  scheme  rather  than  a  standard  programmed  burn  for  the  explosive  detonation. 
Figure  It)  presents  a  comparison  of  the  finite  element  and  experimental  plate  dent  test  results. 

CONCLUSION 

An  advanced  thermodynamic  equation  of  state  applicable  to  problems  involving  overdriven 
detonation  has  been  developed.  The  equation  of  state  has  been  parameterized  for  octol  75/25  and  a 
comparison  of  experimental  results  to  computational  results  has  been  completed.  The  new  equation  of  state 
maintains  required  low  pressure  expansion  behavior  while  providing  a  better  high  pressure  description 
applicable  to  overdriven  detonation.  The  new  equation  of  state  will  be  used  for  evaluating  advanced  multiple 
mode  warheads  concepts  based  on  charge  designs  employing  explosives  with  different  properties  and  initiation 
logic.  Finally,  the  new  equation  of  stale  will  be  used  for  warheads  development  and  analysis  on  several  olhcr 
programs  including  Insensitive  Munitions,  More  Powerful  Explosives  and  Target  Defeat.  Typical  applications 
include  shaped  charge  wave  shaping,  peripheral  initiated  explosively  formed  projectiles,  and  Mach  stem 
detonations. 
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Figure  1.  Pressure  versus  relative  specific  volume  for  the  principle  isentrope  of  octol  75/25  below  (be 
Chapman- Jouguct  state.  The  thermochemical  calculations  (BKWR  and  JCZ3)  agree  fairly  well  with  the 
standard  JWL.  _  _  _  _  .  t 

xio 0  EOS  Calculations 


high  Pressure  Behavior  (Above  CJ) 

Figure  2.  Pressure  versus  relative  specific  volume  for  the  principle  isentrope  and  reactive  Hugoniot  of  octol 
75/25  above  the  Chapman- Jouguct  state.  The  standard  JWL  underpredicts  the  high  pressure  region. 
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Figure  3.  The  Gruneisen  parameter  versus  relative  specific  volume  for  the  principle  isentrope  and  reactive 
Hugoniot. 
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Low  Pressure  Behavior  (Below  CJ) 


Figure  4.  Pressure  versus  relative  specific  volume  for  the  principle  isentrope  of  octol  75/25  below  the 
Chapman- Jouguet  Mate.  The  BKWR  and  JWLB  calculations  agree  fairly  well. 
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Figure  S.  Plots  of  pressure  versus  relative  specific  volume  for  the  principle  isentrope  and  reactive  HugonkX  of 
octol  75/25  above  the  Chapman-Jouguet  state.  The  BKWR  and  JWLB  calculations  agree  very  closely. 


Figure  6.  Plot  of  the  Gruncisen  parameter  versus  specific  volume  for  the  principle  isentrope  and  reactive 
Hugoniot.  The  BKWR  and  JWLB  calculations  agree  very  clcaely. 
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Figure  &  Plate  dent  test  finite  element  computation  initial  and  final  configurations. 
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Figure  9.  Dynamic  finite  element  and  experimental  displacement  versus  time  for  the  cylinder  expansion  test. 
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Figure  10.  Dynamic  finite  element  and  experimental  dent  plate  profile  comparison. 

9 


Cylinder  Test 


DISTRIBUTION  LIST 


Commander 

U.  S.  Army  Armament  Research,  Development  and  Engineering  Center 
A'lTN:  SMCAR-IMI-1  (5) 

SMCAR-CO 
SMCAR-TD 
SMCAR-TDC 
SMCAR-AE  (3) 

SMCAR-AEE  (3) 

Picul  inny  Arsenal,  NJ  07806-5000 

(Commander 

U.  S.  Army  Armament  Research,  Development  and  Engineering  Center 
ATTN:  AMSMC-GCL(D) 

AMSMC-DSM-B(D) 

Picatinny  Arsenal,  NJ  07806-5000 

Commander 

U.  S.  Army  Materiel  Command 
ATrN:  AMCCN-C 

AMCLD.B.  Foneroff 
5001  Eisenhower  Avenue 
Alexandria,  VA  22333-0001 

Commander 
U.  S.  Army  LABCOM 
AIT-N:  TD  R.  Vitali 
2800  Powdermill  Rd 
Adclphi,  MD  20783-1145 

Commander 

U.  S.  Army  Armament,  Munitions  and  Chemical  Command 
A1TN:  AMSMC-DCG 
AMSMC-DSM-A 

Aberdeen  Proving  Ground, MD  21010-5423 
Commander 

U.  S.  Army  Depot  Systems  Command 
ATI’N:  AMSUSUV 
AMSDS-T 
AMSDS-SM-SPA 
AMSDS-CG 

Chambcrsburg.PA  17201-4170 
Administrator 

Defense  Technical  Information  Center 
ATI'N:  Accessions  Division  (12) 

Qimcron  Station 
Alexandria,  VA  22304-6145 

Director 

U.  S.  Army  Materiel  Systems  Analysis  Activity 


11 


ATTN:  AMXSY-MP 
AMXSY-RW 

Aberdeen  Proving  Ground, MD  21005-5066 
Commander 

Chemical  Research, Development  and  Engineering  Center 
U.  S.  Army  Armament,  Munitions  and  Chemical  Command 
ATTN:  SMCCR-MSI-I 
SMCCR-CO 
SMCCR-MU 
SMCCR-SPS-M 
SMCAR-RSP-A 

Aberdeen  Proving  Ground,  MD  21010-5423 
Director 

Ballistic  Research  Laboratory 
ATTN:  AMXBR-OD-ST 
SLCBR-TB-EE 

Aberdeen  Proving  Ground.MD  21005-5066 
Chief 

Benet  Weapons  Laboratory,  CCAC 
Armament  Research, Development  and  Engineering  Center 
U.  S.  Army  Armament,  Munitions  and  Chemical  Command 
ATTN:  SMCAR-CCB-TL 
Watervlict,  NY  12189-5000 

Commander 

U.  S.  Army  Armament, Munitions  and  Chemical  Command 
ATTN:  SMCAR-FSP-D 
SMCAR-DS 
AMSMC-IMF-L 
Rock  Island,  IL  61299-6000 

Director 

U.  S.  Army  TRADOC  Systems  Analysis  Activity 

ATTN:  ATAA-SL 

White  Sands  Missile  Range,  NM  88002 

Commander 
Naval  Weapons  Center 
ATTN:  L  Smith,  Code  320 
A  Amster,  Code  385 
R.  Reed,  Jr.,  Code  388 
K.  J.  Graham,  Code  3835 
China  Lake,  CA  93555 

Director 

Lawrence  Livermore  National  Laboratory 
University  of  California 
P.O.  Box  808 
ATTN:  M.  Finger 
C  Tarver 
L  Green 


12 


Livermore,  CA  94550 


Director 

Ii»  Alamos  National  Laboratory 
ATTN:  L  Hull 

T.  Bennion 
J.  Chapyak 
J.  Ritchie 

Los  Amnios,  NM  87545 
Director 

Sandia  National  Laboratory 
Albuquerque,  NM  87185 

Coi.,mandcr 

Naval  Surface  Weapons  Center 
A'lTN:  Code  G13 
Dahlgrccn,  VA  22448 

Commander 

Naval  Surface  Weapons  Center 
ATFN:  L  Roslund,  R122 
M.  Stosz,  R  121 
Code  X2U,  Library 
E.  Zimct,  R13 
R.R.  Bcrreckcr,  K13 
J.W.  Forbes,  R13 
SJ.  Jacobs,  RIO 
C.  Dickinson 
J.  Short,  R12 
Silver  Spring, MD  20910 

Air  Force  Armament  Laboratory 
ATTN:  J.  Forster,  AFATL-MNW 
G.  Parson,  AFATL-MNE 
Eglin  AFB.FL  32542-5434 

Dyna  East  Corporation 
ATTN:  William  Flis 
3201  Arch  Street 
Philadelphia,  PA  19104 

Aerojet  Ordnance  Company 
A'lTN:  Ditn  Donnli 
2521  Michelle  l)r 
IXtslin,  CA  92680 


